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ABSTRACT 


An  antenna  engineering  handbook  and  database  program  has  been  constructed  by 
engineers  at  the  Lawrence  Livermore  National  Laboratory  (LLNL)  using  the  Numerical 
Electromagnetics  Code  (NEC)  antenna  modeling  program  to  prepare  data  performance 
on  tactical  field  communication  antennas  used  by  the  Army.  It  is  desirable  to  have  this 
information  installed  on  a  personnel  computer  (PC),  using  relational  database  tech¬ 
niques  to  select  antennas  based  on  performance  criteria. 

This  thesis  obtains  and  analyses  current  distributions  and  radiation  pattern  data  by 
using  NEC  for  the  following  set  of  four  (4)  high  frequency  (HF)  tactical  generic  antennas 
to  be  used  in  future  spherical  mode  expansion  work:  a  quarter  wavelength  basic  whip, 
a  one-wavelength  horizontal  quad  Loop,  a  564-foot  longwire,  and  a  sloping  "vee  beam" 
dipole. 

The  results  of  this  study  show  that  the  basic  whip  antenna  provides  good  ground 
wave  communication,  but  it  has  poor  near  vertical  incident  skywave  (N'YIS)  perform¬ 
ance.  The  current  distribution  has  the  characteristics  of  standing  waves.  The  horizontal 
quad  loop  antenna  is  good  for  NV1S  and  medium  range  skywave  communications.  The 
current  distribution  is  sinusoidal  and  continuous  around  the  loop.  The  long  wire  an¬ 
tenna  allows  short,  medium  and  long  range  communications  and  a  standing  wave  cur¬ 
rent  distribution  occurs  along  the  antenna  axis  due  to  non-termination.  The  sloping  "vee 
beam"  antenna  favors  long  range  communication  and  the  current  distribution  is  mainly 
that  of  travelling  sinusoidal  waves.  V 

Because  of  their  well-known  efficiency,  the  basic  whip  and  quad  loop  can  be  used 
as  reference  standards  for  the  spherical  mode  expansion  work.  The  longwire  and  sloping 
"vee  beam"  antenna  are  unwieldy,  but  they  are  effective  as  base  station  antennas. 
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1.  INTRODUCTION 


A.  BACKGROUND 

A  seven  (7)  volume  antenna  engineering  handbook,  has  been  produced  at  the 
Lawrence  Livermore  National  Laborotory  (LLNL)  using  the  Numerical 
Electromagnetics  Code  (NEC)  antenna  modeling  program  for  the  US  Army  Information 
Systems  Engineering  Center  (USAISEC),  Ft.  Huachuca,  Arizona.  [Refs.  1.2]. 

A  host  of  field  tactical  antenna  types  and  configurations  have  been  characterized  for 
optimizing  performance  for  high  frequency  (HF)  communications  when  interfaced  with 
existing  propagation  codes  and  models. 

It  is  desirable  to  have  all  useful  antenna  information  installed  in  a  database  program 
on  a  personnel  computer  (PC)  for  selecting  of  best  candidate  antennas  for  a  given  op¬ 
erational  scencrio.  Initial  work  was  started  on  this  project  several  years  ago  and  has 
continued  to  date  [Ref.  3]. 

Initial  efforts  here  concentrated  on  applying  existing  software  and  new  programs 
necessary  for  completing  this  effort  and  on  producing  a  computer  program  for  the  PC. 

Most  of  the  antennas  in  the  antenna  handbook  are  less  than  a  few  wavelengths  in 
extent,  limiting  the  directivity.  This  limitation  translates  into  patterns  which  are 
dcscribable  in  terms  of  a  few  low  order  spherical  modes. 

A  proposal,  made  by  Dr.  T.  L.  Simpson,  Electrical  and  Computer  Engineering  De¬ 
partment,  University  of  South  Carolina,  was  to  solve  the  problem  of  pattern  generation 
by  finding  a  more  natural  and  simple  way  of  expressing  the  currents  in  terms  of  a 
spherical  modal  expansion.  "Since  the  general  solution  of  the  electromagnetic  wave 
equation  in  spherical  coordinates  can  be  conveniently  expressed  as  a  sum  over  a  set  of 
orthogonal  eigenfunctions,  or  modes  [Refs.  4,  5,  6],  one  can  always  express  a  physically 
realizable  pattern  f  unction  such  as  an  expansion.  In  effect,  the  use  of  the  spatial  filtering 
properties  of  the  antenna  reduces  the  antenna  pattern  to  its  fundamentally  simplest 
form."  [Ref.  7). 

It  is  desired  to  investigate  and  complete  the  representation  of  LIE  antenna  patterns 
using  Spherical  Mode  Expansion.  The  result  of  this  work  will  be  used  to  generate  pat¬ 
terns  on  the  PC  with  minimal  disk  storage  requirements  and  fast  execution  speed. 
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B.  SCOPE  OF  THE  THESIS 


In  order  to  represent  HF  antenna  patterns  using  Spherical  Mode  Expansion,  NEC 
results  for  radiation  pattern  calculations  and  current  distribution  data  are  required  for 
the  support  of  research  in  the  spherical  modal  solutions. 

This  thesis  presents  current  distribution  and  radiation  pattern  data  produced  by 
NEC  for  the  HF  range  of  5  and  30  MHz  for  the  following  set  of  four  (4)  tactical  generic 
antennas:  a  quarter  wavelength  basic  whip  with  8'  ground  stake,  a  one-wavelength  hor¬ 
izontal  quad  loop,  a  564-foot  longwire,  and  a  sloping  "vee  beam"  dipole. 

C.  SELECTION  AND  MODELING  CRITERIA  FOR  THE  ANTENNAS 
1.  Ground  and  Frequency  Selection 

Several  options  are  available  for  modeling  antennas  which  are  located  near  the 
ground.  Either  perfect  ground  or  imperfect  ground  can  be  specified,  with  the  NEC  user 
supplying  values  for  the  conductivity  (<t)  and  relative  permittivity  (£,),  in  the  case  of  finite 
ground,  as  shown  by  the  following  table: 


Table  1.  GROUND  TYPES 


Conductivity 
(mlt os  mi 

Relative 

Permittivity 

Land  Type 

0.00022 

2.5 

flat  desert,  cities 

0.0012 

7.0 

mountains:  steep  rocky  lulls 

0.O03 

10.0 

average  ground 

o.on 

1 3.0 

pastoral  land,  medium  hills 

0.065 

22.5 

rich  farm  land 

0.15 

34.0 

rice  paddy 

The  ground  constants  (conductivity  and  relative  permittivity)  used  throughout 
this  thesis  are  "average"  ground  values  characteristics  of  many  rural  sites,  so  that  the 
performance  calculated  by  the  computer  would  be  similar  to  the  results  one  could  expect 
to  find  under  typical  tactical  conditions. 

The  frequencies  selected  for  this  study  are  5  MHz,  which  is  the  logarithmic 
mean  for  frequencies  over  which  near  vertical  incident  skywave  (NVIS)  communications 
are  conducted,  and  30  MHz,  which  is  the  highest  frequency  defined  by  the  International 

*) 


Radio  Consultative  Commitee  (CCIR)  for  the  HF  band  (3  to  30  MHz).  Ionospheric 
skywave  propagation  is  dominant  at  these  frequencies  and  is  used  for  skyvvave  links. 

The  reason  for  selecting  these  frequencies  is  to  calculate  the  performance  of  the 
generic  antennas  at  the  low  and  high  frequency  limits  of  the  HF  band.  Thus  the  results 
bracket  the  performance  over  the  entire  HF  band. 

2.  Antenna  Selection 

Several  criteria  were  used  when  selecting  the  antennas  for  use  in  the  Antenna 
Handbook.  The  following  performance  criteria  were  selected  as  desirable  for  field  con¬ 
ditions: 

1.  Communication  effectiveness. 

2.  Ruggedness. 

3.  Easy  of  construction  in  the  field. 

4.  Simplicity. 

The  Antenna  Engineering  Handbook  presently  contains  the  following  set  of 
antenna  model  configurations: 

1.  Basic  whip  with  and  without  a  radio  set. 

2.  Basic  half  square. 

3.  Inverted  half  square. 

4.  Folded  half  square. 

5.  3-element  half  square. 

6.  Inverted  3-elemcnt  half  square. 

7.  Horizontal  vee  beam. 

8.  Sloping  vee  beam. 

9.  Sloping  vee  beam,  legs  terminated. 

10.  Horizontal  dipole. 

11.  Inverted- vee  dipole. 

12.  Horizontal  square  loop  (quad). 

13.  Vertical  triangular  loop  (delta  Loop). 

14.  Basic  longwire. 

15.  Sloping  longwire. 

16.  Inverted-vee  longwire. 
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From  these,  a  set  of  four  generic  antennas  mentioned  in  part  B  has  been  se¬ 


lected. 


Besides  military  use,  some  of  the  antennas  selected  for  this  work  are  commonly 
used  by  amateur,  commercial,  and  maritime  HF  stations.  These  were  selected  because 
some  of  them  have  buried  parts,  others  do  not  touch  the  ground,  or  the  physical  sizes 
differ  substantially.  Some  of  them  are  slanted  while  others  have  axes  parallel  to 
cartesian  coordinates,  etc.  Thus  the  simulation  program  can  be  tested  under  different 
conditions  for  these  well-known  antennas. 

D.  DESCRIPTION  OF  THE  ANTENNAS 

1.  Basic  Whip 

•  This  antenna  is  a  quarter  wavelength  long  at  the  operating  frequency.  It  is  verti¬ 
cally  polarized  and  is  fed  against  ground. 

•  The  radius  of  the  antenna  is  0.015  m. 

•  It  has  also  an  S'  ground  stake  buried  beneath  the  ground  as  seen  in  Figure  1. 

2.  One  Wavelength  Horizontal  Quad  Loop 

•  The  One  wavelength  quad  loop  antenna  has  one-quarter  wavelength  sides  and  is 
placed  horizontally  a  quarter-wavelength  high  over  the  ground  plane.  The  feed 
point  is  at  one  corner. 

•  It  has  0.001  m.  radius.  The  structure  is  shown  in  Figure  2. 

3.  564- Foot  Longwire 

•  This  is  a  564-foot  long  wire  antenna  with  a  45°  slanted  feeder,  and  is  placed  40  feet 
high  above  the  ground.  The  feeder  is  driven  by  a  voltage  source  whose  lower  end 
is  tied  to  a  2-foot  ground  stake. 

•  The  radius  of  the  wire  is  also  0.001  m.  The  structure  is  shown  in  Figure  3. 

4.  Sloping  "Vee  Beam"  Dipole 

•  This  antenna  may  be  visualized  as  two  diverging  lines  terminated  through  a 
600-ohm  resistor  to  a  ground  stake.  It  is  fed  at  the  apex. 

•  The  apex  height  is  40  feet  above  the  ground  and  both  ends  are  7  feet  high.  The 
length  of  each  leg  is  52S  feet. 

•  The  conductor  is  also  0.001  m.  in  radius.  The  structure  is  shown  in  Figure  4. 

E.  INTRODUCTION  TO  THE  NUMERICAL  ELECTROMAGNETICS  CODE 

Computer  modeling  of  antennas  started  in  the  late  1960's.  Since  then  it  has  become 
a  powerful  tool  for  antenna  design.  Computer  codes  have  been  developed  based  on  the 
Mcthod-of-Moments,  Geometrical  Theory  of  Diffraction,  or  integration  of  Maxwell's 
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QUARTER  WAVE  WHIP 


THETA  =  60.00  PHI  =  60.00  ETA  =  90.00 


Figure  1.  Quarter  Wavelength  Basic  Whip  with  8'  Ground  Stake 
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Figure  2.  One  Wavelength  Horizontal  Quad  Loop 
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-FOOT  LONGWIRE  ANTENNA  WITH  SLANTED  FEEDER 


.  A 

> 


THETA  =  60.00  PH!  =  60.00  ETA  =  90.00 


Figure  3.  564- Foot  Longwire 
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Equations.  The  Numerical  Electromagnetics  code  -  Method  of  Moments  [Ref.  2]  has 
become  the  most  widely  used  of  all  thin  wire  codes  for  modeling  resonant-sized  anten¬ 
nas.  The  reasons  are  that  NEC  features  systematic  updating  and  extension  of  its  capa¬ 
bilities,  user-oriented  documentation  and  accessibility  of  its  developers  for  user 
assistance. 

NEC  has  been  under  development  for  more  than  ten  years  (earlier  known  as 
BRACT  and  AMP).  It  is  a  hybrid  code  which  uses  an  Electric  Field  Integral  Equation 
(EFIE)  to  model  wire  like  objects  and  a  Magnetic  Field  Integral  Equation  (MFIE)  to 
model  closed  surfaces  with  time  harmonic  excitation.  NEC  is  commonly  applied  to 
modeling  antennas  in  VLF-to-VHF  applications  on  ships,  vehicles  or  on  the  ground. 
A  number  of  features  are  included  for  efficient  modeling  of  antennas  and  scatterers  in 
their  environments,  including  antennas  which  are  interacting  with  or  buried  in  a  finitely 
conducting  ground.  Recently,  a  capability  for  modeling  insulated  wires  in  the  air  or 
ground  has  added  to  NEC  and  improved  accuracy  for  electrically  small  antennas  and  the 
treatment  of  an  abrupt  change  in  wire  radius  [Ref.  8]. 

The  utility  of  NEC  has  been  enchanced  by  the  development  of  an  interactive 
graphics  utility  (IGUANA)  [Ref.  9)  for  developing  the  model  description  and  examining 
the  results.  IGUANA  supports  both  NEC  and  the  Mcthod-of  Moments  Code 
MINTNEC  [Ref.  10]  for  use  on  personnel  computers. 
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II.  PERFORMANCE  RESULTS  OF  ANTENNAS 


A.  GENERAL 

•  All  the  antennas  have  been  excited  by  a  voltage  source  (applied  E  field  source). 

•  The  antennas  are  driven  by  1  kW  input  power. 

•  The  radiation  pattern  range  is  1  km. 

•  The  wavelengths  in  the  average  ground  have  been  calculated  and  found  as  15.6  m 
for  5  MHz,  and  3.1  m  for  30  MHz.  These  data  have  been  used  to  set  segment 
lengths  for  buried  ground  stakes. 

•  All  of  the  calculations  have  been  done  in  double  precision  at  32  bits  for  accuracy. 


B.  BASIC  WHIP  WITH  8'  GROUND  STAKE 

This  antenna  is  basically  a  simple  monopole  with  an  8'  ground  stake.  For  spot  fre¬ 
quency  operation,  a  monopole  antenna  is  satisfactory  at  a  coastal  site,  but  poor  at  low 
conductivity  sites;  its  low  angle  coverage  is  limited  by  the  amount  of  ground  improve¬ 
ment.  Therefore  a  ground  stake  has  been  added  to  the  basic  monopole  structure. 

At  5  MHz,  the  antenna  structure  is  15  m  above  the  ground  (  0.25  ).  ),  and  is  modeled 
by  6  segments,  such  that  all  the  segment  lengths  are  less  than  6  m  (  0.1  /.  ).  For  the 
ground  stake  3  segments  were  used  due  to  the  smaller  wavelength  in  dirt. 

At  30  MHz,  the  length  of  the  antenna  itself  is  2.5  m,  and  it  has  5  segments.  Also 
for  the  buried  part.  15  segments  wer  used  to  model  the  ground  stake. 

At  5  MHz,  where  the  ground  stake  is  electrically  short,  a  part  of  the  total  base  cur¬ 
rent  comes  from  the  return  current  which  flows  in  the  ground  beyond  the  ground  stake. 
The  current  distribution  is  shown  in  Figure  5.  It  is  a  standing  wave  current  distribution. 
The  current  in  the  ground  stake  decreases  rapidly  and  flows  outward  in  the  ground.  Also 
the  phase  is  constant  which  is  a  characteristic  of  standing  waves. 

For  this  current  distribution,  a  radiation  pattern  has  been  calculated  and  plotted  in 
Figure  6.  The  maximum  field  strength  occurs  at  a  27°  take-off  angle. 

The  radiation  resistance  was  calculated  as 


R-rad 


rad 


500 


I  r 


1(2.9  +/0.1)2 


=  58.6  Q. 


(1) 
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At  30  MHz  the  ground  stake  is  about  a  quarter  wavelength  long.  It  provides  good 
ground  improvement  favoring  ground  wave  communication.  The  standing  wave  current 
distribution  is  shown  in  Figure  6.  The  phase  difference  between  the  whip  antenna  and 
the  ground  stake  is  180°  as  expected. 

For  this  frequency,  a  radiation  pattern  has  been  calculated  and  is  shown  in  Figure 
8.  The  maximum  field  strength  occurs  at  29°  take-ofF  angle. 

The  radiation  resistance  at  30  VI  Hz  was  found  as 


Prad  _  500 

|/2|  |(2.9  +/0.1)2| 


50.6  n. 


(2) 


C.  ONE  WAVELENGTH  HORIZONTAL  QUAD  LOOP 

The  one  wavelength  quad  loop  antenna  can  be  considered  also  as  a  modified  folded 
dipole.  By  pulling  the  dipole  wires  apart  at  the  center,  the  single-turn  quad  loop  antenna 
is  obtained.  The  length  of  each  side  is  0.25  . 

At  5  MHz,  the  antenna  is  placed  15  m  above  the  ground.  Each  side  was  modeled 
by  15  segments.  At  30  MHz  the  antenna  is  2.5  m  above  the  ground  and  each  side  was 
modeled  by  13  segments. 

For  a  one  wavelength  perimeter  it  is  reasonable  that  the  current  distribution  is 
sinusoidal  as  shown  in  Figure  9  for  5  MHz  and  in  Figure  13  for  30  MHz.  It  is  contin¬ 
uous  around  the  loop  and  the  magnitude  is  minimum  at  the  corners  between  the  feed 
point  sides  and  the  far  sides.  The  phase  reverses  every  0.5  /  . 

Principal  patterns  at  5  MHz  are  shown  in  the  following  order;  Figure  10  shows  the 
YZ-plane  pattern  plot  of  Ee.  In  this  plane  the  half-power  beamwidth  (HPBW)  is  76°; 
Figure  1 1  shows  the  XZ-plane  pattern  plot  of  Eg  and  the  HPBW  is  1 18°;  Figure  12  shows 
the  XY-plane  pattern  plot  of  £c. 

The  radiation  resistance  was  calculated  as 


Prad 


1000 


1(1.9  +j\.l){ -1.9-71.7)1 


=  142.9  Q. 


(3) 


At  30  MHz  the  YZ-plane  pattern  plot  of  Ee  is  shown  in  Figure  14  with  HPBW  = 
74°;  Figure  15  shows  the  XZ-plane  pattern  plot  of  Ee  with  HPBW  =  116°;  Figure  16 
shows  the  XY-plane  pattern  plot  of  Et.  The  radiation  resistance  is 


R  _ _ 1000 _ 

rad  Kl.9+yi.9X  -1.9  -jl.9)\ 


=  132.3  a 


(4) 
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Current  Phase  Current  Magnitude 


Figure  5 


QUARTER  WAVE  WHIP  OVER  AVE.  GROUND  /  5  MHZ 

WITH  8”  GROUND  STAKE/INCLUDING  THE  GND  STAKE  CURRENT 


Electrical  Distance  (Degrees) 


Current  Distribution  of  Basic  Whip  Antenna  at  5  MHz.:  Bold  line  in¬ 
dicates  the  current  distribution  on  the  ground  stake. 
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QUARTER  WAVE  WHIP  WITH  8’  GROUND  STAKE  OVER  AVERAGE  GND. 
ELEVATION  PATTERN  /  FREQ=5  MHZ. 
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Figure  7.  Current  Distribution  of  Basic  Whip  Antenna  at  30  MHz.:  Bold  line  in¬ 
dicates  the  current  distribution  on  the  ground  stake. 
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QUARTER  WAVE  WHIP  WITH  8'  GROUND  STAKE  OVER  AVERAGE  GROUND. 
ELEVATION  PATTERN  /  FRE0=30  MHZ. 


Figure  8.  Elevation  Pattern  of  Basic  Whip  Antenna  at  30  MHz.: 

c,  =  10.0,  a  —  0.003 
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Current  Phase  Current  Magnitude 


ONE  WAVELENGTH  HORIZONTAL  QUAD  LOOP  /  FREQ=  5  MHZ. 

OVER  AVERAGE  GROUND  /  HEIGHT=  0.25  WAVELENGTH. 


Figure  10.  The  YZ-pIane  pattern  of  Quad  Loop  Antenna  at  5  MHz. 

c,  =  10.0,  a  =  0.003 


ONE  WAVELENGTH  HORIZONTAL  QUAD  LOOP  /  HEIGHT=.25  WAVELENGTH 
OVER  AVERAGE  GROUND  /  ELEVATION  PATTERN  PHI=0  /  FREQ=5  MHZ. 


Figure  11.  The  XZ-p!ane  pattern  of  Quad  Loop  Antenna  at  5  MHz.: 

e,  =  10.0,  o  =  0.003 
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ONE  WAVELENGTH  HORIZONTAL  QUAD  LOOP  /  HE1GHT=.25  WAVELENGTH 
OVER  AVERAGE  GROUND  /  AZIMUTH  PATTERN  /  FREQ=5  MHZ. 


o 


Figure  12.  The  XY-plane  pattern  of  Quad  Loop  Antenna  at  5  MHz.: 


Figure  13.  Current  Distribution  of  Quad  Loop  Antenna  at  30  MHz. 
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ONE  WAVELENGTH  HORIZONTAL  QUAD  LOOP  /  HEIGHT=.25  WAVELENGTH 


OVER  AVERAGE  GROUND  /  ELEVATION  PATTERN  PHI=90  /  FREQ=30  MHZ 


PATTERN  CAIN  IN  081 

...-...■horizontal 


ELEVATION  ANGLE 


Figure  14.  The  YZ-plane  pattern  of  Quad  Loop  Antenna  at  30  MHz. 

Cr  =  10.0.  0  =  0.003 
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ONE  WAVELENGTH  HORIZONTAL  QUAD  LOOP  /  HEtGHT=.25  WAVELENGTH 
OVER  AVERAGE  GROUND  /  ELEVATION  PATTERN  PHI=0  /  FREQ=30  MHZ. 


Figure  IS.  The  XZ-plane  pattern  of  Quad  Loop  Antenna  at  30  MHz.: 

c,  =  10.0,  a  =  0.003 
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ONE  WAVELENGTH  HORIZONTAL  QUAD  LOOP  /  HEIGHT=.25  WAVELENGTH 
OVER  AVERAGE  GROUND  /  AZIMUTH  PATTERN  /  FREQ=30  MHZ. 


o 


Figure  16.  The  XY-plane  pattern  of  Quad  Loop  Antenna  at  30  MHz.: 

z,  =  10.0,(7  =  0.003 
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D.  LONG-WIRE  ANTENNA 

The  wire  antennas  introduced  so  far  have  been  resonant-size  antennas.  The  wave 
travelling  outward  and  reflected  from  the  end  sets  up  a  standing-wave  type  of  current 
distribution.  From  this  section  on,  non-resonant  size  antennas  will  be  introduced. 

A  long-wire  antenna  is  one  that  is  greater  than  one  halfwave  length  long.  When 
such  antennas  are  unterminated  at  the  end,  the  current  distribution  is  mainly  that  of  a 
standing  wave.  If  the  antenna  is  cut  to  a  resonant  length,  the  input  impedance  will  be 
resistive.  The  theoretical  current  distribution  is  a  sinusoidal  standing  wave.  But  because 
of  the  loss  due  to  radiation,  the  actual  radiation  pattern  differs  from  the  theoretical  cal¬ 
culation.  The  first  recognizable  difference  is  that  the  lobes  tilt  toward  the  unfed  end. 
This  difference  between  actual  patterns  and  theoretical  patterns  (based  on  standing 
waves  only)  is  much  less  in  the  case  of  center-fed  antennas. 

If  a  longwire  is  terminated  in  a  resistance  equal  to  its  characteristic  impedance,  the 
current  distribution  is  mainly  that  of  a  travelling  wave.  The  longer  the  wire,  the  smaller 
will  be  the  angle  between  the  wire  and  the  main  lobe.  The  important  difference  between 
the  patterns  of  terminated  and  unterminated  wire  antennas  is  the  absence  of  large  rear 
lobes  in  the  terminated  case. 

In  this  study  the  long  wire  was  564  feet  long  and  placed  40  feet  above  the  ground. 
It  was  unterminated. 

At  5  VlIIz  the  longwire  is  on  the  order  of  3  wavelengths  (171.9  m)  and  was  modeled 
by  51  segments.  Also  this  antenna  has  a  45°  slanted  feeder  segmented  by  8  pieces.  This 
feeder  was  driven  by  a  voltage  source  whose  lower  end  is  tied  to  a  2-foot  ground  stake. 
This  ground  stake  was  radius  tapered  and  had  2  segments.  At  30  MHz  the  antenna  is 
approximately  17  wavelengths  long  and  has  200  segments.  The  slanted  feeder  was 
modeled  by  45  segments  and  the  ground  stake  had  4  segments. 

For  the  unterminated  case  with  a  pure  sinusoidal  current  distribution  (standing 
wave),  the  source  current  can  be  expressed  as  a  superposition  of  two  uniform  waves 
travelling  in  opposite  directions  as 

/(. x)  =  /0  sin[m7r(  -j-  +  0.5)]  (5) 

Lk  ■  • 

where  m  is  an  integer  equal  to  — — ,  L  is  the  length  of  the  antenna,  and  k  is  the 

~>jr 

wavenumber  (-^4-). 

/ 
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Figure  17  shows  the  current  distribution  at  5  MHz.  Standing  waves  occur  along  the 
antenna  axis  due  to  non-termination.  Also  the  magnitude  decreases  slightly  due  to  ra¬ 
diation.  The  phase  also  reverses  every  half  wavelength. 

Figure  18  shows  the  elevation  pattern  at  5  MHz.  It  is  symmetrical  about  the  Z-axis 
and  has  a  main  beam  calculated  by 


a 


m 


cos  '(1 


0.371 

m  ’ 


(6) 


which  is  30°,  and  additional  beams  in  even.'  a  =  180  —  o.m  direction  providing  good  NVIS, 
medium  range,  and  long  range  communications. 

Figure  19  shows  the  azimuth  pattern  for  a  30°  take-off  angle.  Figure  20  shows  the 
same  pattern  for  a  60°  take-off  angle.  Figure  21  shows  the  azimuth  pattern  for  an  80° 
take-off  angle.  All  the  radiation  patterns  above  are  symmetric  along  the  antenna  axes 
and  the  vertical  gains  are  dominant  to  the  horizontal  gains. 

Radiation  resistance  at  this  frequency  was  calculated  as 


Rrad~ 


500 


I(0.S+/1.1)‘| 


247.4  Q. 


(7) 


Figure  22  shows  the  current  distribution  at  30  MHz.  Standing  waves  exist  as  in  the 
former  case  but  the  effect  of  losses  is  higher  than  for  the  same  antenna  at  5  MHz.  In 
addition,  there  are  also  a  few  travelling  waves.  The  minimum  magnitude  of  the  current 
occurs  in  the  middle  of  the  antenna. 

Figure  23  shows  the  elevation  pattern  at  30  MHz.  It  has  a  main  beam  at  13°.  Al¬ 
though  the  radiation  pattern  has  symmetry  about  the  z-axis,  the  gain  of  the  front  lobe 
is  larger  than  the  back  lobe. 

Figure  24  and  25  show  the  azimuth  patterns  for  15°  and  25°  take-off  angles  respec¬ 
tively.  Both  radiation  patterns  show  the  main  lobe  in  the  forward  direction  larger  than 
in  the  reverse  direction.  This  indicates  also  the  existence  of  travelling  waves  on  the  an¬ 
tenna. 

Radiation  resistance  at  this  frequency  was  calculated  as 

Rrad= - -°°  —=171.1  Q.  (8) 

1(1.1  -j\-2Y\ 
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564-FOOT  LONGWIRE  WITH  SLANTED  FEEDER  /  FREQ=  5  MHZ. 
OVER  AVERAGE  GROUND  /  HEIGHT  =  40  FEET. 
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Figure  18.  Elevation  Pattern  of  Long  Wire  Antenna  at  5  MHz.: 

z,  =  10.0.  a  =  0.003 
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Figure  19.  Azimuth  Pattern  of  Long  Wire  Antenna  at  5  MHz.  a  =  30°: 
t,  —  1 0.0.  a  -  0.003 
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564-FOOT  LONGWiRE  WITH  SLANTED  FEEDER  /  HEIGHT=  40  FEET. 
OVER  AVERAGE  GROUND  /  AZIMUTH  PATTERN  THETA=30  /  FREQ=5  MHZ. 
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564-FOOT  LONGWIRE  WITH  SLANTED  FEEDER  /  HEIGHT=4Q  FEET. 
OVER  AVERAGE  GROUND  /  AZIMUTH  PATTERN  THETA=10  /  FREQ=5  MHZ. 


o 


Figure  21.  Azimuth  Pattern  of  Long  Wire  Antenna  at  5  MHz.  a  =  80°: 
cr  =  10.0,  a  =  0.003 
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Current  Phase  Current  Magnitude 


Figure  22.  Current  Distribution  of  Long  Wire  Antenna  at  30  MHz. 


31 


Figure  23.  Elevation  Pattern  of  Long  Wire  Antenna  at  30  MHz.: 


c,  =  10.0,  a  =  0.003 
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Figure  24.  Azimuth  Pattern  of  Long  Wire  Antenna  at  30  MHz.  a=15c 

c,  =  10.0,  a  =  0.003 
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564-FOOT  LONGWIRE  WITH  SLANTED  FEEDER  /  HEIGHT=40  FEET. 
OVER  AVERAGE  GROUND  /  AZIMUTH  PATTERN  THETA=65  /  FREQ=30  MHZ 


o 


Figure  25.  Azimuth  Pattern  of  Long  Wire  Antenna  at  30  MHz.  a  =  25°: 
t,  =  10.0,  o  =  0.003 
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E.  SLOPING  VEE  BEAM  ANTENNA 

This  is  a  travelling-wave,  balanced  antenna.  The  angle  between  the  descending  wires 
was  chosen  so  that  the  principal  lobes  of  two  radiators  are  superimposed.  The  mast 
height  is  chosen  to  make  the  antenna  tilt  angle  the  same  as  the  take-off  angle.  Half  of 
the  radiated  power  is  reflected  from  the  ground.  Therefore  the  ground  should  be  smooth 
and  free  from  obstructions.  This  antenna  is  horizontally  polarized  and  has  a  relatively 
high  input  impedance. 

The  modeled  antenna  had  two  528-foot  descending  wires.  The  interior  angle  be¬ 
tween  two  wires  was  60°.  The  feed  point  was  40  feet  above  the  ground  and  each  leg  of 
the  "Vee"  was  terminated  through  a  600  Q  resistor  to  a  2  feet  ground  stake. 

At  5  MHz  each  528-foot  leg  has  45  segments.  The  load  is  placed  on  one  segment 
at  this  frequency.  The  ground  stake  was  modeled  by  2  segments. 

At  30  MHz  each  leg  had  221  segments.  The  wire  connecting  the  load  to  the  ground 
stake  has  4  tapered-length  segments.  Each  ground  stake  was  modeled  with  4  segments. 
The  gain  at  lower  take-off  angles  provides  good  long  range  communication.  Figure  33 
shows  the  azimuth  pattern  at  a  10°  take-off  angle.  Although  it  is  horizontally  polarized, 
some  of  the  radiated  power  comes  from  ground  reflection. 

For  the  current  distribution  on  antenna  wire  No.  1, 

I,  -  (9) 

and  for  the  current  distribution  on  antenna  wire  No.  2, 


where  In  is  the  current  amplitude  and  s  is  the  distance  from  the  feed  point. 

Figure  26  shows  that  the  current  distribution  at  5  MHz  is  essentially  a  travelling 
wave  sinusoidal  current  distribution.  The  magnitude  decreases  gradually  toward  the 
load  and  the  phase  reverses  even  half  wavelength. 

Figure  27  shows  the  elevation  pattern  of  Ee  at  <f>  =  60°,  the  lobe  created  by  one 
descending  wire.  Figure  28  shows  the  YZ-plane  pattern  plot  of  Ee.  The  lobe  nearest  to 
the  antenna  wire  is  the  largest  (the  main  lobe)  and  all  others  progressively  diminish  in 
strength.  Figure  29  shows  the  azimuth  plane  pattern,  Et,  at  an  elevation  angle  of  20°  . 
The  horizontal  pattern  gain  in  the  direction  of  antenna  is  maximum  and  the  back  lobes 
are  introduced  by  ground  reflections.  The  input  impedance  is  highly  resistive.  At  this 
frequency  the  radiation  resistance  was  calculated  as 
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(10) 


p  _  A  VW  _ 

rad~  |(1.1  H-J0.006X  — 1.1  — >0.006)1  ” 

Figure  30  demonstrates  that  the  current  distribution  at  30  MHz  is  also  a  sinusoidal 
distribution  with  logarithmically  decreasing  magnitude. 

At  30  MHz  the  antenna  is  about  16  wavelengths  long.  Figure  31  shows  the  ele¬ 
vation  pattern  at  4>  =  60°.  Figure  32  shows  also  the  YZ-plane  pattern  of  Et.  Here  the 
principles  lobes  are  superimposed.  The  gain  is  maximum  at  a  low  take-off  angle  pro¬ 
viding  good  long  range  communication.  Figure  33  shows  the  azimuth  pattern  at  a  10° 
take-ofF  angle. 

At  this  frequency  radiation  resistance  was  calculated  as 


R-rad 


1000 


1(1.2  +/0.005)(  -1.2  -j0.005)| 


=  690  a 


(11) 
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Figure  27.  Elevation  Pattern  of  Sloping  Vee  Antenna  at  5  MHz.  <t>  =  60° 

c,  —  10.0,  <7  =  0.003 
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SLOPING  "VEE  BEAM"  DIPOLE  /  LENGTH  OF  EACH  LEG  -  528  FEET 
OVER  AVERAGE  GROUND  /  ELEVATION  PATTERN  PHI=90  /  FREQ=5  MHZ. 


90 


180 


PATTERN  GAIN  IN  OBI 

.  .........horizontal 

•  —  VERTICAL 

- . .  TOTAL 

elevation  angle 


Figure  28.  YZ-Plane  Pattern  of  Sloping  Vee  Antenna  at  5  MHz. 

£,  =  10.0,  <r  =  0.003 


SLOPING  "VEE  BEAM"  DIPOLE  /  LENGTH  OF  EACH  LEG  =  528  FEET 
OVER  AVERAGE  GROUND  /  AZIMUTH  PATTERN  THETA=70  /  FREQ=5  MHZ. 


Figure  29.  Azimuth  Pattern  of  Sloping  Vee  Antenna  at  5  MHz.  a  =  20°: 
£,  =  10.0,  <7  =  0.003 
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Current  Phase  Current  Magnitude 


SLOPING  "VEE  BEAM”  DIPOLE  /  ONE  LEG  /  FREQ=3 

OVER  AVERAGE  GROUND  /  INCLUDING  THE  GND.  STAKE  CUF 


Figure  30.  Current  Distribution  of  Sloping  Vee  Antenna  at  30  MHz. 


Figure  31.  Elevation  Pattern  of  Sloping  Vee  Antenna  at  30  MHz.  <t>  =  60c 

c,=  10.0,  a  =  0.003 
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SLOPING  "VEE"  BEAM  DIPOLE  /  LENGTH  OF  EACH  LEG  =  528  FEET. 
OVER  AVERAGE  GROUND  /  ELEVATION  PATTERN  PHI=90  /  FREQ=30  MHZ 


Figure  32.  YZ-Plane  Pattern  of  Sloping  Yee  Antenna  at  30  MHz. 

e,  -  10.0,  a  =  0.003 


Figure  33.  Azimuth  Pattern  of  Sloping  Vee  Antenna  at  30  MHz.  a  =  10° 

e,  =  10.0,  a  =  0.003 
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III.  CONCLUSIONS  AND  RECOMMENDATIONS 
A.  CONCLUSIONS 

This  thesis  investigated  the  performance  of  4  tactical  HF  field  antennas  at  5  and  30 
MHz  over  average  ground  conditions. 

The  basic  whip  antenna  provides  good  ground  wave  communication.  The  gain  of 
the  ground  wave  is  proportional  to  the  ground  conductivity  improvement.  The  skywave 
is  not  as  dependant  on  ground  constants.  It  has  poor  NVIS  performance,  but  its  low 
take-off  angle  favors  long  distance  communication.  At  HF  the  usefulness  of  the  basic 
whip  is  limited  by  the  wide  operating  bandwidth  usually  required.  The  current  distrib¬ 
ution  has  the  characteristics  of  standing  waves. 

The  horizontal  quad  loop  antenna  has  highest  gain  at  high  take-off  angles.  There¬ 
fore  it  has  very  poor  ground  wave  performance.  The  radiation  patterns  shew  that  it  is 
good  for  WIS  and  medium  range  skywave  communications.  The  one  wavelength  quad 
loop  has  the  same  elevation  and  azimuth  pattern  as  a  half  wavelength  dipole,  but  it  has 
more  gain.  Also  the  loop  antenna  has  lower  terminal  resistance  than  a  folded  dipole. 
The  current  distribution  is  sinusoidal  and  continuous  around  the  loop. 

The  long  wire  antenna  allows  short,  medium,  and  long  range  communications.  The 
radiation  patterns  at  different  frequencies  show  that  as  the  frequency  increases,  the  main 
lobe  becomes  more  lobey  and  has  better  vertical  gain  at  lower  angles,  providing  long- 
range  communication.  The  performance  of  the  antenna  remains  practically  unchanged 
within  10  MHz  of  the  design  center  frequency.  Generally,  antennas  which  use  termi¬ 
nated  lines  appear  to  have  more  radiation  resistance  than  unterminated  ones.  Radiation 
and  wire  losses  and  the  presence  of  ground  and  obstacles  alfect  the  pattern  of  any  ideal 
source.  If  improvement  is  required,  it  is  better  to  lengthen  the  longwire  than  to  raise  it. 
The  height  of  the  antenna  can  be  adjusted  so  that  the  reflected  wire  combines  in  phase 
with  the  direct  wave  in  the  direction  of  maximum  radiation.  The  important  consider¬ 
ation  here  is  the  phase  shift  at  the  point  of  reflection  [Ref.  11].  Since  it  is  unterminated, 
a  standing  wave  current  distribution  occurs  along  the  antenna  axis. 

For  the  sloping  "vee  beam"  antenna,  it  is  also  true  that,  like  a  longwire,  this  antenna 
has  also  a  relatively  wide  bandwidth.  Since  half  of  the  radiated  power  reflects  from 
ground,  the  ground  should  be  free  from  obstacles.  The  level  of  the  sidelobes  change 
considerably,  depending  on  the  type  of  ground.  The  level  of  the  largest  sidelobes  is  only 
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about  4  -  7  dB  below  the  main  beam.  This  is  why  the  usefulness  of  the  antenna  is  re¬ 
stricted.  Its  low  take-off  angle  provides  very  good  long-range  communication.  The 
current  distribution  is  mainly  that  of  travelling  sinusoidal  waves. 

As  a  result,  because  of  their  well-known  efficiency,  the  basic  whip  and  quad  loop 
can  be  used  as  reference  standards  for  the  spherical  mode  expansion  work.  The  longwire 
and  528-foot  sloping  "vee  beam  "  dipole  are  large,  unwieldy  and  difficult  to  construct, 
but  they  are  effective  as  base  station  antennas. 

B.  RECOMMENDATIONS 

The  configurations  of  different  antennas  can  be  studied  under  different  conditions 
or  for  different  frequencies  to  provide  more  information  for  spherical  modal  expansion 
work,  if  needed. 
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APPENDIX  A.  NEC  AND  AUXILIARY  PROGRAMS 


A  brief  description  of  some  of  the  programs  found  in  the  NEC  FORTRAN  library 
are  described  below: 

1.  SOMNTX  FOR  SOMMERFELD  /  NORTON  GROUND  METHOD 

This  option  can  be  requested  on  the  GN  card.  This  method  is  available  for  wires 
only.  It  uses  an  exact  solution  for  the  fields  in  the  presence  of  ground  and  is  accurate 
close  to  and  beneath  the  ground.  NEC  reads  interpolation  tables  from  the  file  TAPE21. 
This  input  file  must  be  created  by  running  SOMNTX.  SOMNTX  reads  4  parameters 
from  a  data  file: 

•  EPR  =  relative  dielectric  constant  of  ground  (  t, ) 

•  S1G  =  conductivity  of  ground  in  mho,  m  (  a  ) 

•  FMHZ  =  frequency  in  MHz. 

•  1PT  =  print  option  of  interpolation  tables. 

2.  THE  NUMERICAL  GREEN  S  FUNCTION  (NGF)  OPTION 

The  main  purpose  of  the  NGF  is  to  avoid  unnecessary  repetition  of  calculations. 
Another  reason  for  using  the  NGF  option  is  to  exploit  partial  symmetry  in  a  structure. 
This  option  can  be  requested  by  a  \VG  card  and  a  GF  card. 

3.  DNPGNEC 

A  double  precison  version  of  NEC-  ?  available  in  mainframe  at  NPS,  it  is  good 
for  modeling  structures  up  to  300  segment.  DN’PGIOOO,  another  version  of  NEC,  is 
capable  of  handling  up  to  1000  segments. 

4.  PLOT  PROGRAMS 

•  PLOTDGLP  —  Plots  data  geometry  to  a  laser  printer. 

•  PLOTDG79  —  Plots  data  geometry  to  a  3279  terminal. 

•  PLOTRPHM  --  Plots  horizontal  and  vertical  components  of  radiation  patterns. 

Use  this  to  get  azimuth  plots  in  the  horizontal  plane. 

•  PLOTRPVM  --  Plots  radiation  pattern  in  the  vertical  plane. 

•  CURRPLOT  —  This  program  plots  current  magnitute  and  phase  for  selected  seg¬ 
ments  and  tags  or  selected  groups  of  segments  or  tags. 
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5.  RVAL 

It  is  used  for  variable  segmentation  or  tapered  radius  or  both.  RVAL  asks  the 
length  of  the  first  and  the  last  segments  and  the  total  length  of  the  wire.  Then  it  calcu¬ 
lates  the  segment  lengths  and  the  segmentation  ratio.  It  is  used  on  GC  cards,  which  are 
read  only  when  the  radius  on  a  GW  card  is  set  to  zero. 

6.  DIRTLAM 

This  program  calculates  the  wavelength  in  dirt.  The  result  is  used  to  set  segment 
lengths  for  wires  buried  under  ground.  It  asks  the  relative  dielectric  constant, 
conductivity  and  operation  frequency  and  gives  the  wavelength  in  the  ground. 
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APPENDIX  B.  NEC  DATA  SETS  FOR  CURRENTS 

CM  GEOMETRY  :  QUARTER  WAVELENGTH  BASIC  WHIP  WITH  8'  GROUND  ROD 

CM  RADIO  :  NONE 

CM  FREQUENCY  :  5  MEGAHERTZ 

CM  GROUND  :  EPSILON  =  10;  SIGMA  =  0.  003 

CE 

GW  1,6,0,0,0,0,0,14.989,0.015 
GW  2,3,0,0,0,0,0,-2.438,0.015 
GE  0,2 

FR  0,0, 0,0,  5 

GN  2,0,0,0,10,. 003 

EX  0,1,1,01,342. 755 

PL  1,2, 0,0 

XQ 

WG 

EN 


CM  GEOMETRY  :  QUARTER  WAVELENGTH  BASIC  WHIP  WITH  8'  GROUND  ROD 

CM  RADIO  :  NONE 

CM  FREQUENCY  :  30  MEGAHERTZ 

CM  GROUND  :  EPSILON  =  10;  SIGMA  =  .003 

CE 

GW  1,5,  0,0,0,  0,0,2.498,  0.015 
GW  2,15,  0,0,0,  0,0,-2.438,  0.015 
GE  0,2 

FR  0,0, 0,0,  30, 

GN  2, 0,0,0,  10,.  003 
PL  1,2, 0,0 
EX  0,1,1,01,336.641 

XQ 

WG 

EN 


CM  GEOMETRY  :  ONE  WAVELENGTH  HORIZONTAL  QUAD  LOOP 

CM  HEIGHT  =  0. 25  WAVELENGTH 

CM  CONDUCTOR  =  #12  AWG 

CM  FREQUENCY  :  5  MEGAHERTZ 

CM  GROUND  :  EPSILON  =  10;  SIGMA  =  .003 

CE 

GW  1,12,  0,0,14.989,  -10.599,10.599,14.989,  0.001 
GW  2,12,  -10.599,10.599,14.989,  0,21.198,14.989,  0.001 
GW  3,12,  0,0,14.989,  10.599,10.599,14.989,  0.001 
GW  4,12,  10.599,10.599,14.989,  0,21.198,14.989,  0.001 
GE  -1,2 

FR  0,0, 0,0,  5, 

GN  2, 0,0,0,  10,. 003 
PL  1,2, 0,0 

EX  0,1,1,01,  251. 115 
EX  0,3,1,01,  -251. 115 
XQ 
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WG 

EN 


CM  GEOMETRY  :  ONE  WAVELENGTH  HORIZONTAL  QUAD  LOOP 

CM  HEIGHT  =  0. 25  WAVELENGTH 

CM  CONDUCTOR  =  # 12  AWG 

CM  FREQUENCY  :  30  MEGAHERTZ 

CM  GROUND  :  EPSILON  =  10;  SIGMA  =  .  003 

CE 

GW  1,13,  0,0,2.498,  -1.766,1.766,2.498,  0.001 
GW  2,13,  -1.766,1.766,2.498,  0,3.  532,2.  498,  0.001 
GW  3,13,  0,0,2.498,  1.766,1.766,2.498,  0.001 
GW  4,13,  1.766,1.766,2.498,  0,3.532,2.498,  0.001 
GE  -1,2 

FR  0,0, 0,0,  30, 

GN  2, 0,0,0,  10,. 003 
PL  1,2, 0,0 

EX  0,1,1,01,  252. 345 

EX  0,3,1,01,  -252.345 

XQ 

WG 

EN 


CM  GEOMETRY  :  564-FOOT  LONGWIRE  WITH  45  DEGREE  SLANTED  FEEDER 
CM  HEIGHT  ABOVE  GROUND  =  40  FEET 

CM  CONDUCTOR  IS  #  12  AWG 

CM  FEEDER  IS  DRIVEN  BY  A  VOLTAGE  SOURCE  WHOSE  LOWER 

CM  END  IS  TIED  TO  A  2 -FOOT  GROUND  STAKE 

CM  FREQUENCY  :  5  MHZ 

CM  GROUND  :  EPSILON  =  10  ;  SIGMA  =  .  003 
CE 

GW  1,2,  0.0762,0.06096,-0.6096,  0.0762,0.06096,0 
GC  0,0,  1,  0. 004,0. 002 

GW  2,1,  0.0762,0.06096,0,  0.0762,0.06096,1.302,  0.001 

GW  3,8,  0.0762,0.06096,1.302,  -10.8138,0.06096,12.192,  0.001 

GW  4,51,  -10.8138,0.06096,12.192,  -182.721,0.06096,12.192,  0.001 

GE  -1,2 

FR  0,0, 0,0,  5 

GN  2, 0,0,0,  10,.  OP* 

PL  1,2, 0,0 

EX  0,3,1,01,  1204. 350 

XQ 

WG 

EN 


CM  GEOMETRY  :  564-FOOT  LONGWIRE  WITH  45  DEGREE  SLANTED  FEEDER 
CM  HEIGHT  ABOVE  GROUND  =  40  FEET 

CM  CONDUCTOR  IS  #  12  AWG 

CM  FEEDER  IS  DRIVEN  BY  A  VOLTAGE  SOURCE  WHOSE  LOWER 

CM  END  IS  TIED  TO  A  2 -FOOT  GROUND  STAKE 

CM  FREQUENCY  :  30  MHZ 

CM  GROUND  :  EPSILON  =  10  ;  SIGMA  =  .  003 
CE 

GW  1,4,  0.0762,0.06096,-0.6096,  0.0762,0.06096,0 
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GC  0,0,  1,  0. 004,0. 002 

GW  2,6,  0.0762,0.06096,0,  0.0762,0.06096,1.302,  0.001 
GW  3,45,  0.0762,0.06096,1.302,  -10.8138,0.06096,12.192,  0.001 
GW  4,200,  -10.8138,0.06096,12.192,  -182.721,0.06096,12.192,  0.001 
GE  -1,2 

FR  0,0, 0,0,  30 
GN  2, 0,0,0,  10,. 003 
PL  1,2, 0,0 

EX  0,3,1,01,  872.360 

XQ 

WG 

EN 


CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CE 

GW 

GW 

GW 

GC 

GW 

GW 

GW 

GC 

GE 

LD 

LD 

FR 

GN 

PL 

EX 

EX 

XQ 

WG 

EN 


GEOMETRY  :  SLOPING  "VEE  BEAM"  DIPOLE 
APEX  HEIGHT  =  40  FEET 
BOTH  ENDS  AT  HEIGHT  =  7  FEET 
LENGTH  OF  EACH  LEG  OF  "VEE”  =528  FEET 
EACH  LEG  OF  "VEE"  IS  TERMINATED  THROUGH  A  600-0 HM 
RESISTOR  TO  A  2 -FOOT  GROUND  STAKE 
INTERIOR  ANGLE  OF  "VEE"  =  60  DEGREES 
CONDUCTOR  IS  #  12  AWG 
FREQUENCY:  5  MHZ 

GROUND  :  EPSILON  =  10  ;  SIGMA  =  . 003 

1.45,  0,0,12.192,  80.3099,139.1008,2.1336,  0.001 

2.1,  80.3099,139.1008,2.1336,  80.3099,139.1008,0,  0.001 

3.2,  80.3099,139.1008,0,  80.3099,139.1008,-0.6096 
0,0,  1,  0.002,0.004 

4.45,  0,0,12.192,  -80.3099,139.1008,2.1336,  0.001 

5.1,  -80.3099,139.1008,2.1336,  -80.3099,139.1008,0,  0.001 

6.2,  -80.3099,139.1008,0,  -80.3099,139.1008,-0.6096 
0,0,  1,  0.002,0.004 

-1,2 

4,  2,1,1,  600 
4,  5,1,1,  600 
0 ,0 , 0 , 0 ,  5 
2, 0,0,0,  10,. 003 
1,2, 0,0 

0,1,1,01,  436.485 
0,4,1,01,  -436.485 


CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CE 


GEOMETRY  :  SLOPING  "VEE  BEAM"  DIPOLE 
APEX  HEIGHT  =  40  FEET 
BOTH  ENDS  AT  HEIGHT  =  7  FEET 
LENGTH  OF  EACH  LEG  OF  "VEE"  =  528  FEET 
EACH  LEG  OF  "VEE"  IS  TERMINATED  THROUGH  A  600 -OHM 
RESISTOR  TO  A  2 -FOOT  GROUND  STAKE 
INTERIOR  ANGLE  OF  "VEE"  =  60  DEGREES 
CONDUCTOR  IS  #  12  AWG 
FREQUENCY:  5  MHZ 

GROUND  :  EPSILON  =  10  ;  SIGMA  =  .003 
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GW  1,221,  0,0,12.192,  80.3099,139.1008,2.1336,  0.001 
GW  2,4,  80.3099,139.1008,2.1336,  80.3099,139.1008,0, 

GC  0,0,1. 306,0. 001,0. 001 

GW  3,4,  80.3099,139.1008,0,  80.3099,139.1008,-0.6096 
GC  0,0,  1,  0. 002,0. 004 

GW  4,221,  0,0,12.192,  -80.3099,139.1008,2.1336,  0.001 
GW  5,4,  -80.3099,139.1008,2.1336,  -80.3099,139.1008,0, 
GC  0,0,1. 306,0. 001,0. 001 

GW  6,4,  -80.3099,139.1008,0,  -80.3099,139.1008,-0.6096 
GC  0,0,  1,  0.002,0.004 
GE  -1,2 

LD  4,  2,4,4,  600 
LD  4,  5,4,4,  600 
FR  0,0, 0,0,  30 
GN  2, 0,0,0,  10,. 003 
PL  1,2, 0,0 

EX  0,1,1,01,  415.428 

EX  0,4,1,01,  -415.428 

XQ 

WG 

EN 
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APPENDIX  C.  NUMERICAL  METHODS  AND  TECHNIQUES  FOR 

ANTENNA  ANALYSIS 

A.  THE  METHOD  OF  SOLUTION 

NEC  uses  both  an  electric-field  integral  equation  (EF1E)  and  a  magnetic-field  inte¬ 
gral  equation  (MFIE).  The  EFIE  in  NEC  is  well-suited  for  thin-wire  structures  of  small 
or  vanishing  conductor  volume.  MFIE  is  only  useful  for  closed  structures,  especially 
those  having  large  smooth  surfaces.  For  a  structure  containing  both  wires  and  surfaces, 
the  EFIE  and  MFIE  are  coupled  in  NEC. 

To  model  wires,  the  thin  wire  approximation  is  used  in  the  EFIE  to  reduce  it  to  a 
scalar  integral  equation.  The  assumptions  used  are  [Ref.  8  ]  : 

1.  Transverse  currents  are  neglected  relative  to  axial  currents  on  the  wire. 

2.  The  circumferential  variation  in  the  axial  current  is  neglected. 

3.  The  current  is  represented  by  a  filament  on  the  wire  axis. 

4.  The  boundary  condition  on  the  electric  field  is  enforced  in  the  axial  direction  only. 

These  approximations  are  valid  for  the  wire  radius  much  less  than  the  wavelength 
and  much  less  than  the  wire  length.  An  alternative  kernel  for  the  EFIE,  based  on  an 
extended  thin-wire  approximation,  is  also  included  in  NEC  for  wires  having  too  large  a 
radius  for  the  thin-wire  approximation  [Ref.  12).  With  the  thin-wire  approximation  the 
EFIE  becomes: 


-  s  .  E'(F)  = 


-J 
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where 

s 

s' 

E'(F) 

co 

/(s') 

£ 


=  unit  vector  along  the  wire  axis. 
=  distance  along  the  wire  axis. 

=  incident  electric  field  at  r. 

=  2,t/,  where  f  is  frequency. 

=  axial  current. 

=  permittivity. 
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k 


=  cojuc  =  phase  constant. 

7  =  observation  point. 

T  =  source  point. 

[Jjkr  i 

— — —  =  free  space  Green's  function. 

R  =  (F-F') 

The  MFIE  for  a  closed  surface  S  is: 


Js(7')  x  Vgds'rzs'  (13) 

where 

1,(7)  =  surface  current  density. 

Hmc(7)  =  incident  magnetic  field  at  the  observation  point. 

h  =  unit  vector  normal  to  the  surface  at  7. 

The  integral  equations  (12)  and  (13)  are  solved  numerically  in  NEC  by  a  form  of  the 
method  of  moments.  This  solution  includes  expanding  the  unknown  current  in  a  sum¬ 
mation  of  basis  functions  with  unknown  amplitudes  and  enforcing  equality  of  weighted 
integrals  of  fields.  The  weighting  functions  for  both  wires  and  surfaces  in  NEC  are 
chosen  to  be  delta  functions  which  resulted  in  a  point  sampling  of  the  fields.  Wires  are 
divided  into  short  segments  with  a  sample  point  at  the  center  of  each  segment.  The  sur¬ 
faces  are  approximated  by  a  set  of  fiat  patches  or  facets  with  a  sample  point  at  the  center 
of  each  patch. 

Delta  functions  are  also  used  in  NEC  as  the  current  expansion  functions  on  a  sur¬ 
faces.  For  the  MFIE,  Galerkin's  method  [Ref.  13]  provides  good  accuracy  on  large 
smooth  surfaces.  Due  to  nature  of  the  integral-equation  kernels,  however,  the  choice 
of  current  expansion  functions  is  more  critical  in  the  EFIE  than  the  MFIE. 

A  more  realistic  representation  of  the  surface  current  is  needed  where  a  wire  con¬ 
nects  to  a  surface  than  the  delta  function  expansion  normally  used  with  the  MFIE.  In 
the  region  of  the  wire  connection,  the  surface  current  includes  a  singular  component  due 
to  the  current  flowing  from  the  wire  onto  the  surface. 

As  a.  result  of  the  continuity  conditions  on  current  and  charge  there  remains  one 
unknown  associated  with  each  basis  function  to  be  determined  by  solving  the 
electromagnetic  interaction  equations.  Enforcing  appropriate  conditions  on  the  current 


Js(7)  =  2n  x  Hlnc(7)  + 


and  charge  results  in  rapid  convergence  of  the  solution.  The  conditions  may  be  inap¬ 
propriate  in  some  cases,  such  as  for  short  segments  at  a  change  in  radius  and  at  junc¬ 
tions  of  closely  spaced  wires.  The  result  is  then  slower  convergence.  More  accurate 
conditions  for  these  cases  are  currently  being  developed  for  NEC  [Ref.  8]. 

B.  CAPABILITIES  OF  NEC 

There  are  basically  two  versions  of  NEC  written  in  FORTRAN.  N'EC-2  has  no 
distribution  restrictions  and  limited  a  maximum  of  300  segments  to  model  a  given 
structure.  On  the  other  hand  NEC-3  is  classified  as  Defense  Critical  Technology  and 
can  be  released  only  to  DoD  agencies  and  their  contractors  without  specific  DoD  ap¬ 
proval.  NEC-3  can  be  used  up  to  300,  1000,  or  2000  segments,  depending  on  the  ver¬ 
sions  of  NEC  being  used.  The  capabilities  of  these  codes  are  summarized  below. 

1.  Input 

A  user-oriented  data  card  set  for  a  single  run  consists  of  three  types  of  data 
cards.  The  deck  must  begin  with  one  or  more  comment  cards  which  can  contain  a  brief 
description  and  structure  parameters  for  the  run.  These  are  followed  by  geometry  data 
cards  which  permits  defining  straight  wires,  arcs,  and  surfaces.  Electrical  connections 
are  determined  in  the  program  by  searching  for  wire  ends  and  patch  centers  that  coin¬ 
cide.  Shifts,  rotation,  and  reflections  can  be  used  in  building  structures.  Finally,  a  sec¬ 
tion  of  program  control  cards  specifies  electrical  parameters  such  as  frequency,  loading 
and  excitation,  and  requests  calculation  of  antenna  currents  and  fields. 

2,  Output 

Outputs  are  selectable  by  input  parameters  and  may  include  one  or  more  fol¬ 
lowing  data: 

1.  charge  density  on  wires. 

2.  current. 

3.  input  impedance,  admittance,  and  power. 

4.  radiated  power,  ohmic  loss,  efficiency. 

5.  radiated  field,  power  gain  and  directive  gain. 

6.  average  gain. 

7.  near  E  and  H  fields. 

8.  maximum  coupling. 

9.  receiving  patterns. 

10.  scattering  cross  sections. 
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3.  Source  modeling 

A  structure  may  be  excited  by  one  of  the  following  sources: 

1.  voltage  source  (applied-E-field  source). 

2.  incident  plane  wave,  linear  polarization. 

3.  incident  plane  wave,  right-hand  elliptic  polarization. 

4.  incident  plane  wave,  left-hand  elliptic  polarization. 

5.  elementary  current  source. 

6.  voltage  source  (current-slope-discontinuity). 

4.  Nonradiating  Networks  and  Transmission  lines 

Nonradiating  two-port  networks  and  transmission  lines  may  connect  points  on 
wires.  The  defining  parameters  are  characteristic  impedance,  length,  and  shunt 
admittance. 

5.  Loading 

Lumped  or  distributed  RLC  loads  may  be  specified  on  wires.  The  options  that 
the  user  can  specify  are: 

1.  series  RLC,  input  ohms,  henries,  farads. 

2.  parellel  RLC,  input  ohms,  henries,  farads. 

3.  series  RLC,  input  ohms  meter,  henries  meter,  farads/meter. 

4.  parallel  RLC.  input  ohms  meter,  henries  meter,  farads  meter. 

Also  the  conductivity  of  a  round  wire  may  be  specified  and  the  impedance 
computed  taking  account  of  skin  depth. 

6.  Ground  Effects 

For  an  antenna  structure  in  or  near  the  ground,  there  are  three  options  avail¬ 
able.  A  perfectly  conducting  ground  is  modeled  by  including  the  image  field  in  the  kernel 
of  the  integral  equations.  An  approximate  model  for  a  finitely  conducting  ground  uses 
the  image  modified  by  Fresnel  reflection  coefficients.  This  approximation  is  usable  for 
antennas  at  least  0.1  to  0.2  wavelengths  above  the  ground. 

NEC-3  includes  an  accurate  treatment  for  wire  structures  above,  below,  or 
penetrating  the  ground  surface.  The  solution  is  based  on  the  Sommerfeld  integral  for¬ 
mulation  for  the  field  near  the  interface. 
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7.  Modeling  Guidelines 

Some  of  the  guidelines  for  modeling  are  summarized  below. 

For  accurate  results,  the  lengths  of  wire  segments  should  be  less  than  about  0.1 
).  .  Longer  segments  up  to  about  0.14  /  may  be  acceptable  on  long  straight  wires  or 
noncritical  parts  of  a  structure.  For  critical  regions,  shorter  segments  0.05  /.  or  less  may¬ 
be  needed.  The  minimum  segment  length  set  by  the  computational  precision  is  about 
10'M.  In  double  precision  the  minimum  length  is  about  10~8/. 

The  wire  radius,  a,  relative  to  is  limited  by  approximations  used  in  the  kernel 
of  the  EFIE.  Two  options  are  available,  the  thin-wire  kernel  and  the  extended  thin-wire 
kernel.  With  the  extended  thin-wire  kernel  the  ratio  of  segment  length  to  radius  can  be 
as  small  as  0.5.  With  the  thin-wire  approximation,  the  limit  is  2.  In  either  case,  only- 
currents  in  the  axial  direction  on  a  segment  are  considered,  and  there  is  no  allowance  for 
variation  of  the  current  around  the  wire  circumference.  The  acceptability  of  these  ap¬ 
proximations  depends  on  both  the  value  of  a  7.  and  the  tendency  of  the  excitation  to 
produce  circumferential  currents.  Discontinuites  in  wire  radius  should  be  kept  to  less 
than  a  factor  of  two. 

Electrically  small  loops  present  special  accuracy  problems.  In  single  precision 
the  minimum  loop  circumference  is  about  0.002  /  ;  while  in  double  precision  it  is  about 

io-4;. . 

The  surfaces  modeled  with  patches  must  be  closed.  A  reasonable  maximum  for 
the  area  of  a  surface  patch  appears  to  be  0.04  square  wavelengths. 

An  important  consideration  in  using  NEC  is  the  solution  time  versus  model  size 
since  this  may  limit  the  amount  of  detail  that  can  be  modeled  and  the  segment  and  patch 
densities.  For  a  model  using  X  wire  segments,  the  solution  time  in  seconds  on  a  CDC 
computer  is  approximately  [Ref.  8  ] 

T  =  3(10-J)A.V2/.U  +  2(K )~6).\3IM2  (14) 

where  M  is  the  number  of  degrees  of  symmetry  and  k  is  1  for  free  space,  2  for  perfectly- 
conducting  ground  modeled  with  the  reflection  coefficient  approximation,  and  4  to  S  for 
Sommerfield  integral  treatment  of  finitely  conducting  ground. 

The  first  ferm  in  Eq.  (3)  is  the  time  to  fill  the  interaction  matrix  and  the  second 
term  is  the  time  to  factor  it  into  triangular  parts  for  solution.  For  a  model  using  Np 
patches,  the  solution  time  is  about  [Ref.  8] 

T=  (Kr5)k(2\p)2IM  +  2(lQ-6)(2\/l.\f2  (15) 
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since  two  rows  and  columns  in  the  matrix  are  associated  with  each  patch. 

When  the  Sommerfeld  ground  treatment  is  used,  a  fixed  time  of  about  15  sec¬ 
onds  [Ref.  8]  is  needed  to  generate  the  interpolation  tables.  The  tables  depend  only  on 
the  ground  parameters  and  frequency,  and  can  be  saved  for  use  in  any  case  in  which 
these  parameters  are  the  same. 

C.  CONCLUSION 

NEC  is  a  versatile  code  for  modeling  antennas  and  their  environment  including 
transmission  lines,  networks,  loading  and  ground  effects.  It  is  written  in  FORTRAN 
and  consists  of  about  10.000  lines  of  code. 

Recently  released  extensions  of  NEC  include  a  model  for  insulated  wires  and  im¬ 
proved  accuracy  in  modeling  electrically  small  antennas  and  discontinuities  in  wire  ra¬ 
dius.  A  long-range  goal  is  to  incorporate  expert-systems  technology  to  simplify  the 
process  of  building  models  that  conform  to  the  modeling  guidelines  and  to  verify  that 
the  code  is  producing  correct  results  [Ref.  S]. 

NEC  will  continue  to  be  a  widely  used  antenna  modeling  program,  due  to  its  con¬ 
venient  operation,  documentation,  availability,  and  continuing  support. 
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